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px in central Au+Au collisions. The thermal-shower component of the recombination 
is found to be important for K and A, but only in a minor way for <p and Vt in the 



intermediate to high px region. The normalization and inverse slope of the thermal 
partons in the strange sector are determined by fitting the low-p^ data. At higher px 



^ , the data of K, A, 4> an d O in the log scale are all well reproduced in our study that 



extends the thermal contribution and includes the shower contribution. The calculated 
result on the A/K ratio rises to a maximum of around 2 at px ~ 4 GeV/c, arching 
over the data in linear scale. The production of <fi and f2 are shown to arise mainly 
from the recombination of thermal partons, thus exhibiting exponential pt dependences 
in agreement with the data. Their ratio, Rq/^, rises linearly to px ~ 4 GeV/c and 
develops a maximum at px ~ 5.5 GeV/c. It is argued that the px spectra of cj> and O 
reveal directly the partonic nature of the thermal source that characterizes quark-gluon 
plasma. Comments are made on the f2 puzzle due to the simultaneous observation of 
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both the exponential behavior of the f2 spectrum in px and the existence of \ow-px 
particles associated with O as trigger. 

1 Introduction 

The production of strange particles has always been a subject of great interest in heavy- 
ion collisions because of their relevance to possible signatures of deconfinement and flavor 
equilibration [Q [2]. Strangeness enhancement that has been observed at various colliding 
energies is a phenomenon associated with soft particles in the bulk matter [31 H]. At high 
transverse momentum (px), on the other hand, the production of jets does not favor strange 
particles, whose fragmentation functions are suppressed compared to those for non-strange 
particles. At intermediate px range between the two extremes the px distribution depends 
sensitively on both the strangeness content and the production mechanism. It has been shown 
that in that px range the spectra of 7r and p can be well described by parton recombination 
[5]. In this paper we study the production of the openly strange particles K, A and Q, as well 
as the hidden strange hadron 0. The aim is to get a broad view of strangeness production 
at intermediate and high px in a unified framework. 

The four particles emphasize four different aspects of hadronization: K and A involve 
non-strange quarks, whereas (f) and Q do not; K and <ft are mesons, whereas A and Q are 
baryons, although <p and A are comparable in masses. To get all four spectra correctly at 
all px would require a high degree of coordination in the theoretical description. We have 
already a fairly good description of the production of hadrons in the non-strange sector [5]; 
however, for the strange sector we must deal with a new set of problems. Gluon conversion 
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feeds the s quarks in the thermal medium and leads to strangeness enhancement (6] . On the 
other hand, the production of s quarks at high pt in jets is suppressed. The interplay between 
the enhanced thermal s partons and the suppressed shower s partons results in very unusual 
outcome. Furthermore, how the supply of s quarks is partitioned into various channels of 
hadrons containing open and hidden strangeness is an aspect of the recombination process 
that was not encountered in the non- strange sector, but was studied algebraically by quark 
counting in Refs. [7J [8] without gluon conversion. Our consideration here treats the problem 
in a much wider scope where the distributions of all four strange particles are investigated. 

We shall use the basic recombination formula to compute the px distributions of all four 
particles. The thermal s quark distribution is determined by fitting the low-pr data of K 
production. The shower parton distributions (SPD) are known [9]. Thus one can proceed to 
the calculation of the A spectrum essentially without adjustable parameters. In the cases of 
and Q the competition for multiple strange quarks to form such high mass states would 
lead to suppression. Moreover, since the s quarks from the showers are suppressed compared 
to non-strange shower partons, cf) and Q produced at large px do not benefit from hard 
scattering as much as K and A do, for which light shower partons can contribute. The 
consequence is that the thermal partons play a dominant role — to the extent that the 
production of and Q provides a window to look into the thermal nature of the quark-gluon 
plasma before hadronization. All these problems will be examined quantitatively in this 
paper. 
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2 Formulation of the Problem 

We shall assume that all hadrons produced in heavy-ion collisions are formed by recombina- 
tion of quarks and/or antiquarks, the original formulation of which is given in [TOl [TT] for pp 
collisions. In recent years the recombination model has been extensively studied by many 
groups [121 D2J Ull E] with great success in reproducing the pt spectra in the intermediate px 
region of Au+Au collisions. In our ID description of the recombination process the invariant 
inclusive distribution of a produced meson with momentum p is 

dN M fdp 1 dp 2 

p —7— = / F q q,(p 1 ,p 2 )R M (p l ,p 2 ,p), (1) 

dp J pi p 2 

and for a produced baryon 

dN B fdp 1 dp 2 dp 3 

P —r- = / F qg , g „(p 1 ,p 2 ,p 3 )R B (p 1 ,p 2 ,p 3 ,p). (2) 

dp J pi P2 P3 

The properties of the medium created by the collisions are imbedded in the joint quark 
distributions F qq i and F qq i q ii. The recombination functions (RFs) Rm and R B depend on the 
hadron structure of the particle produced, but for a high mass hadron the normalization can 
in addition be influenced by the competition among possible channels. In the valon model 
description of hadron structure [TTJ [15], the RFs are 

R M (pi,P2,p) = 9m ViV2Gm{vuV2): (3) 
Rb(pi,P2,P3,p) = 9b yi y2yzG B (yi, 2/2,2/3), (4) 

where yi = Pi/p, and gu and g B are statistical factors. Gm and G B are the non-invariant 
probability densities of finding the valons with momentum fractions yi in a meson and a 
baryon, respectively. Equations ([1]) and ([2]) apply for the produced hadron having any 
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momentum p. We consider p only in the transverse plane, and write pt as p so that dN/pxdpT 
becomes (pp )" 1 times the right-hand sides of Eqs. (CD) and ([2j). 

The wave function G(yi) in terms of the momentum fractions of the valons cannot be 
probed for K, <fr, A, and Q, as one can for proton. However, we can make reasonable estimates. 
For K and A, there is a mixture of strange and non-strange constituents with different masses, 
so the corresponding momentum fractions are different. The wave function of K has been 
considered before [TB] ; it has the form 

= W( ■ I kTT\ iftMvi + V2 - 1). (5) 
B[a + 1,0 + 1) 

The ratio of the average momentum fractions y% and y 2 of the u and s type valons, respec- 
tively, should be roughly the ratio of their masses mjj /ms ~ 2/3. That puts a constraint on 
the parameters a and 6, which is 6 = (3a + l)/2. Phenomenological analysis of the process 
K + + p — > 7r + + X in the valon model then results in a = 1 and 6 = 2 [16]. For the hyperon 
A the wave function can be written in the form [B] 

G A (yu y 2 , y 3 ) = B{a + l ^ a + ^ 2)B{a + l ^ +l) 2/1 * (yi + V> + Vs ~ 1) , (6) 

where y\ and yi refer to the non-strange valon, and y 3 to the strange valon. The ratio of 
the average momentum fractions yi/yz is (a + 1)/(/3 + 1), which yields a similar constraint: 
/3 = (3a + l)/2. We suppose that their values are roughly the same as for proton [15], so 
they shall be set at a = 1 and (3 = 2. 

For cf) and Q only one type of valon is involved, so their valon distributions have the 
simple symmetric forms 

G 4> (y l ,y 2 ) = [B(c+l 7 c + l)]- 1 (y l y 2 ) c 8(y l + y 2 -l), (7) 
G n {y u y 2> y 3 ) = [B( 1 + l, 1 + l)B( 1 + l,2 1 + 2)]-\y 1 y 2 y 3 y5(y 1 +y 2 + y 3 -l). (8) 
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Since the masses of those two particles are nearly at the thresholds of the sums of the 
constituent quark masses, loose binding means narrow wave functions in the momentum 
space. We shall consider two cases: (a) c = 7 = 5, and (b) c = 7 = 4. For case (a) or larger 
c and 7, the distributions in Eqs. ([7j) and (jHJ) are narrow enough so that they are not too far 
from the ^-function approximations, i.e., 

G <j> (y 1 ,y 2 ) = 5( Vl - 1/2) 5(y 2 - 1/2), (9) 
G n (3/i,2/ 2> y 3 ) = %i - 1/3) % 2 - 1/3) 5(|/3 -1/3). (10) 

To calculate the pt distributions, we start with the production of K®, for which we 
consider the ds and sd contributions to F qq -i and write 

F dsM = T% + %{S d + St) 12 + TS S + {SS S } , (11) 

where a sum over all relevant quark species is implied when not indicated explicitly. The 
thermal and shower parton distributions of the light quarks are, respectively, [5] 

dN th 
dpi 

S(p 2 )=£E/ dkkfiWSifa/k). (13) 

i Jk o 

The parameters C and T have been determined in [5]: 

C = 23.2 GeV" 1 , T = 0.317 GeV . (14) 

Si is the SPD for a light quark in a shower initiated by a hard parton i, fi(k) is the transverse- 
momentum (k) distribution of hard parton % at midrapidity in central Au+Au collisions, and 
£ is the average fraction of hard partons that can emerge from the dense medium to hadronize. 
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As in [5], fi(k) is taken from Ref. [IT], k is set at 3 GeV/c, and £ is found to be 0.07. For 
the strange quarks the corresponding thermal and shower parton distributions are 

dN th 

%{Pi) = Pi—r~ = C s p 1 exp(-p 1 /T s ), (15) 
dpi 

S S (P2)=CT, dkkf t (k)S°(p 2 /k), (16) 

i Jk o 

where C s and T s are now new parameters. The SPDs are known from [5] and cannot be 
altered. Indeed, Si and S* are so well determined from the meson fragmentation functions 
(FFs) that we have recently been able to calculate the baryon FFs from them in agreement 
with data [TS] , thus solidifying the validity of our approach to fragmentation in the framework 
of recombination. The thermal s quark distribution in Eq. (|T5l) is determined by fitting the 
low-pr data on production. Since K* decays totally into irK, we have g K = 1. 
The last term of Eq. f[TT|) symbolizes 

{SS S } = £Y, dkkfi(k) 

i Jk o 

where {■••,••■} on the RHS represents symmetrization of pi and p^- The two terms in the 
square brackets are not the same because of the valence quark contribution to the (^-initiated 
jet. Despite the appearance of the product form on the LHS, there is only one integral over 
k because both shower partons belong to the same jet inititated by one hard parton %. Since 
we sum i over u, d, s, u, d, s and g, Eq. (fTTl) contains 28 terms. However, when they are 
combined with the RF, they become the FF Di'{p/k) 0. 
For A production we need 

F uds = TTT S + TT S S + TTS S + T S {SS} + T{SS S } + {SSS S } , (18) 

where the last term involves the symmetrization of all three shower partons in the same jet. 
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Clearly, there are now many more terms for F u d s . We shall find that up to px ~ 4 GeV/c, 
the relative importance of each of the six terms in Eq. (TI8]) diminishes in the order they 
appear on the RHS. Nevertheless, they all have to be calculated in order to determine the 
emergence of certain terms, as is increased. For the RF of A the statistical factor is 
1/4. With the parameter a in the A wave function expressed in Eq. ([6]) set at 1 to simulate 
the proton, there are no adjustable parameters in this A production problem. 

For and Q the pt distributions can be written out explicitly if the naively simple forms 
in Eqs. (ED and (TTOT) for their RFs are used. They are 



pdp ppo 



F 8S (p/2,p/2), (19) 



^ = ^F sss (p/3,p/3,p/3), (20) 
pap ppo 

where 

F sS = T S T S + T S S S + {S S S S } } (21) 
Fsss — T S T S T S + T S T S S S + T S {S S S S } + {5 s 5 s iS s } . (22) 

In our calculation we shall, however, use the RFs given in Eqs. (J7J) and (jSJ) and consider 
the two cases given below those equations. We allow the normalization factors g$ and ga 
to be adjustable in order to account for the medium effect on the recombination of ss and 
sss. Because of the limited supply of s quarks and the abundance of the light quarks in 
the vicinity, it is expected that the formation of the multi-strange particles at higher masses 
would be suppressed. Without a principle from outside the recombination model, there is 
no way to determine the suppression factors. In the following we vary and gn as free 
parameters to fit the normalizations of the corresponding spectra. Current data on the 
spectra of strange particles extend to px as high as 6 GeV/c [18j-[25]. 



3 Results 

We first show the result of our calculation for K® production. In Fig. 1 is shown the contri- 
butions from the three components indicated in the figure. The thermal-shower component 
includes both the T S S and TS S subcomponents, although the latter subcomponent is much 
smaller because the shower s quark is suppressed. The contribution from T S S is more im- 
portant because thermal s quarks are not negligible, as evidenced by the TT S component 
that is dominant for px < 3 GeV/c. The {iS5 s } component does not become important until 
p T > 8 GeV/c. 

In fitting the low-p^ region for p? < 2 GeV/c we use 

C s = 13.9 GeV" 1 , T s = 0.33 GeV . (23) 

This value of T s is slightly higher than T = 0.317 GeV for light quarks. Since the low 
Pt phenomenology includes the flow effect, which we do not separately put in, that small 
difference in the inverse slopes is sensible in view of the mass difference involved. Note that 
it is because of the importance of the contribution of the thermal-shower recombination for 
Pt > 3 GeV/c that the px distribution bends up from the exponential behavior to agree 
with the data [23] . 

Having determined C s and T s , we can now proceed to the calculation of A production 
without adjusting any parameters. There are, however, two cases to consider that are related 
to the valence quarks in the showers. The question is whether the leading parton is likely to 
form a leading meson in a jet or a baryon with higher mass, remembering that there is only 
one valence quark in a quark-initiated shower. We calculate the px distributions of A for 
the two cases, where (a) only the sea quarks in the showers are used for the formation of A, 

9 



and (b) all contributions including all valence quarks are included. The results are shown in 
Fig. 2(a) and (b), respectively. There is a small but perceptible difference between the two 
cases for px > 4 GeV/c. The inclusion of valence quarks in (b) naturally has the effect of 
increasing the yield at high px- 

In Fig. 2 we first note that the recombination of three thermal partons (TTT) is dominant 
up to about 4 GeV/c, where the TTS contribution becomes more important at higher px- 
The TSS component is unimportant below 6 GeV/c, and SSS is even less important. Fig. 
2(a) gives a more satisfactory reproduction of the data [23] than Fig. 2(b). It implies that 
the valence quarks that have larger momentum fractions in a jet are mainly in the mode of 
formation of mesons instead of hyperons. Gluon jets make the dominant contribution to the 
shower partons that form the A. 

With the calculated px distributions of K® and A at hand, we can take their ratio and 
compare it with the data on A/K® [221 E3J. Since the latter is shown in linear scale, it is 
far more sensitive to the details of the px distributions than the spectra themselves plotted 
in log scale. What is shown by the solid line in Fig. 3 is the result on the A/K® ratio using 
the px distributions shown in Figs. 1 and 2(a). Our calculated result in Fig. 3 does not 
reproduce the data at large px, although they agree well for p T < 2 GeV/c. The solid line 
has the right structure of a peak reaching almost as high as 2. The disagreement with data 
is an amplification of the small discrepancies between theory and experiment in Figs. 1 and 
2(a), where our calculated distribution for K® in the region 4 < px < 6 GeV/c is slightly 
lower than the data, but for A it is slightly higher than the data. That intermediate px 
region is where the shower contributions become important. Evidently, our model for the 
production of such strange particles is not good enough to yield better than 30% accuracy. 
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When the shower parton contribution is not strong, it is necessary to rely on the accuracy 
of the thermal parton distribution extrapolated to very high pr, higher than can perhaps be 
justified with the simple exponential formula used. Moreover, it is likely that hadronization 
into different strange hadronic channels may not be independent, so the partition problem 
in the strange sector may have to be taken into account. 

For the production of and O we have done the calculation for both cases (a) c = 7 = 5 
and (b) c = 7 = 4 for the valon distributions given in Eqs. ([7]) and (jHJ) that are used in the 
RFs in Eqs. Q and (pEj). For clarity we show only the results for case (a) in Figs. 4 and 5, 
since the other case will be shown to be ruled out by the 0/0 ratio. If we ignore all but the 
first terms on the RHS of Eqs. (12~T1) and (j22j) . we would have only the thermal contributions, 
which from 1115)) . (1201) and (jI51) givtt 

dN * Cl^e-^ , (24) 



pdp s Ap 

dN n ^3 P 
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where C s and T s are as given in Eq. ff23l) . The results from using c = 7 = 5 in Eqs. ([7]) and 
(jHJ) are similar, but cannot be expressed in simple analytic form as in Eqs. ff24|) and (125]) ; 
they are shown by the dashed lines in Figs. 4 and 5. Being from the thermal contributions 
only, they are nearly exponential for px > 2 GeV/c. The other terms in Eqs. (I2ip and 
(I22p involving shower partons give rise to contributions shown by the other lines in those 
figures, which are lower. Their sums represented by the solid lines differ only slightly from 
the dashed lines. It is clear then from Figs. 4 and 5 that the pt distributions of both 



1 These simple formulas are the result of using the <5-function approximations of the wave functions in 
Eqs. (J9]) and (fTO]) for the recombination functions. 
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and Q are dominated by the recombination of thermal partons (T s ) in the regions where 
data are available [2T L 124 1 [25]. The contribution of thermal-shower recombination (T S S S ) to 
(f> production does not exceed that from T S T S until pr > 6 GeV/c. For Q the cross-over is 
above 8 GeV/c, not shown in Fig. 5. The production of shower s quark is suppressed, but 
it is enough to make a few percentage contribution at the high end of pt where data exist c 
The overall normalizations of the distributions are fitted by adjusting the values of g<$> and 
gn to be 

94> = 0.34 , g n = 0.012 . (26) 

They represent the degree of suppression of the rate of recombination of an s quark with an 
s (to form 0) and two other s quarks (to form Q) in the environment of light quarks that 
present competing channels for the formation of kaons and hyperons. The phenomenon is 
similar to that in the A production where the leading s quark prefers to form K rather than 
A, as shown in the two parts of Fig. 2. 

To calculate the ratio Rn/</> we can start with the simplest by using Eqs. ( J24l) and (1251) 
and get 



R%M - (27) 



which is linearly rising. However, that is for thermal-parton recombination with 5-function 
RFs. For ) RFs with c = 7 = 5, the result for thermal contribution only is shown 

by the dashed line in Fig. 6, which is the ratio of the dashed lines in Figs. 4 and 5, and is 
approximately linearly rising. When all terms are included, the ratio of the solid lines in 



2 This result differs from the earlier one in [26] . where a number of approximations were made, and a 
smaller values of £ was used by mistake. 
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Figs. 4 and 5 results in the solid line in Fig. 6. There is a definitive bending down of the 
ratio for pt > 5 GeV/c, although the solid lines in Fig. 4 and 5 deviate only slightly from 
the dashed lines in the log scale. The linear scale in which the ratio is shown magnifies the 
deviation arising from the contributions of the shower partons. The solid line exhibits the 
trend of the data that departs from linearity at pt ~ 4 GeV/c [25]; however, our model 
calculation gives a peak at pt ~ 5 GeV/c which is slightly higher than the data, and with 
the maximum at Rq/c/, ~ 0.25. The decrease of tt/<f> at high p? is due to the increase of 
the spectrum that is lifted up from the exponential behavior by the contribution from 
thermal-shower recombination. That upward bending of the <p spectrum occurs earlier in p? 
than that of Q, with the consequence that tt/<f> decreases for px > 5 GeV/c. In Fig. 6 we 
also show in dashed-dotted line the Q/4> ratio for c = 7 = 4 in the wave function of <fi and 
Q [see Eqs. flTJ) and (|SJ)], corresponding to broader widths in the momentum-fraction space. 
The result is far from the data and should be discarded from further consideration. It does 
confirm our earlier statement that narrow width in the momentum space is expected due to 
the weak binding of constituent quarks in and Q. 

4 Discussion 

In the approach to the hadronization problem that we have adopted there are certain prop- 
erties that we cannot change, while there are others that we can determine only by fitting 
the data. The former involve the hard partons and their shower parton distributions. Of 
course, they can be improved by going to higher order calculations, which we rule out for 
the present consideration. The latter involve the soft thermal partons, for which several 
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parameters have been used for non-strange and strange quarks. Since our aim has never 
been to model the soft component, it does not matter how many parameters have been used 
to reproduce the low pr behavior. Our aim has been to understand strange particle produc- 
tion at intermediate and high p?, and we have found some valuable information to provide 
significant insight into the strange sector. The recombination model has been effective in 
illuminating the physics involved. 

From K® and A spectra we have seen how the deviation from the exponential behavior 
that describes the thermal components can be understood as being due to the emergence of 
the shower contribution arising from jets. It happens at around px ~ 4 GeV/c. Those are 
strange hadrons containing light quarks, which are easier to produce at high pt than the 
strange quarks that are suppressed. The thermal source provides the s quarks for recombi- 
nation. The A/K ratio has a maximum at p? ~ 4 GeV/c, a bit higher than the location of 
the peak in the data [23]. The cause of the decrease of the ratio is the relative ease of K 
over A in acquiring a boost in yield by thermal-shower recombination at high p?. Quanti- 
tatively, we do not expect better than 30% accuracy in the spectra themselves, since lowest 
order pQCD was used in the calculation of hard scattering, and the virtuality of the shower 
partons is not changed according to the p? of the fragmenting hard partons. Despite that 
limitation the fact that the spectra can be reproduced over such a wide range of yield for p T 
up to 6 GeV/c (a span of 5 orders of magnitude) is sufficient to be regarded as support for 
the validity of the recombination model used. 

In the case of <fi and Q the absence of light quarks in them deprives them of the light 
shower partons at high px, so their productions are dominated by the thermal strange-quark 
recombination. The contribution from hard parton scattering is negligible unless pt is large. 
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The small upward bend of the and Q spectra above straight-line behaviors in the log plots 
(i.e., exponential) results in the f2/0 ratio to bend down above pt ~ 5 GeV/c. Q is harder 
to be lifted up by shower s quark than because there are two thermal s quarks instead of 
one, given that shower s partons are suppressed in the first place. 

The significance of the f2/0 ratio in Fig. 6 is not so much the bending down of the solid 
line at high as the rising of the dashed line at intermediate pt in agreement with the data. 
That rising behavior can be seen clearly in the simple analytic expressions in Eqs. (I24I) and 
(I2"5"j) . from which one gets the ratio to be proportional to p (which is abbreviation for p? in 
Fig. 6). Although those expressions are simplified approximations, they exhibit two features: 
(a) the extra factor of p in Eq. ( |25l) compared to ( |24l) is due to the difference in the RFs of 
the baryon versus meson shown in Eqs. (T4|) and ([31), respectively, and (b) the more important 
aspect of those two expressions is the common exponential factors that are cancelled in the 
ratio — those factors are dynamical in origin, while the prefactors are kinematical. The fact 
that the dashed line in Fig. 6 well describes the data up to pt ~ 4 GeV/c implies that both 
and Q have the same dynamical origin of thermal partons, which are assumed to have 
exponential behavior given in Eq. ffTBT) . The surprise is the validity of that exponential form 
out to p ~ 2 GeV/c for observed at pt ~ 4 GeV/c, remembering that T s is only 0.33 
GeV. Thus we have a very robust system of thermal partons consisting of not only the light 
quarks but also the strange quarks. That is the main characteristic of quark-gluon plasma 
[1, 2J. It takes the production of and Q to reveal cleanly this thermal source because their 
exponential behavior at intermediate px is not covered up by the contamination of shower 
partons from hard scattering, as is the case with other mesons and baryons. In that sense 
and Q serve as windows through which one can observe directly the properties of the partonic 
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source. 

The discussion above concerns the shapes of the px distributions of <ft an d fi. The 
normalizations have been adjusted to fit the data by giving the parameters g$ and ga the 
values shown in Eq. fl26|) . The fact that they are less than 1 implies suppression. We interpret 
that as being due to the competition for s quarks in the environment of light quarks during 
hadronization. Evidently, it is less likely for an s to find an s to form than a u or d to form 
K; it is even more difficult to find two other s quarks to form Q. To study the constraints 
arising from competing channels is outside of the scope of this paper, since our emphasis has 
been on the px dependence of the strange hadron production. The values of and gn can, 
however, be useful at a later point when the competiton problem is considered. 

The phenomenological fact that the pt distribution of Q is essentially exponential up to 
p T ps 6 GeV/c leads one to conclude, even without the benefit of our model calculation, 
that Q is produced from the thermal source. It is then natural to predict that any particle 
produced in association with Q must also be from the same thermal source, and hence indis- 
tinguishable from the background. Such a prediction was, in fact, made in the unpublished 
version of this work which stimulated an experimental search for the partners of Q 
above background. The result of that study was reported at Quark Matter 2006, indicating 
that unidentified charged particles were found on the near side of the A0 distribution in 
association with f2 as trigger for 2.5 < p% lg < 4.5 GeV/c and 1.5 < p^ soc < p^ lg [27]. Those 
are rather low values of p^ g and p^ ssoc , lower than what can unambiguously give rise to 
the usual jet structure. For p^ lg < 4.5 GeV/c we see in Fig. 5 that the contribution from 
thermal-thermal-shower recombination is minuscule compared to that of thermal partons 
only. Thus there are no apparent jets to account for the associated particles. This dilemma, 
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called the Q puzzle, was addressed at Quark Matter 2006 and resolved conceptually by the 
notion of phantom jet [28J. Without entering into the details of the possible resolution, we 
mention here only the essence of the idea. It is conjectured that both the fl and its associated 
particles are produced from the ridge induced by jets, the nature of which has been studied 
for other triggers [27] . The ridge is an enhanced thermal source, and is itself not jet-like, but 
without a jet (created by gluon or light quark) there would not be any ridge [29]. A ridge 
without peak may therefore be regarded as the shadow of a phantom jet. If the associated 
particles are in the ridge, then they can be observed above the background. Quantitative 
implementation of this idea remains to be carried out, but probably not until more data 
become available on ridge characteristics. At present, without At] distribution there is no 
identification of a ridge due to the low statistics in Q production. 

In conclusion, we remark that intermediate pt physics in the strange sector at RHIC 
has been substantially illuminated by this study where hadron production up to pr ~ 6 
GeV/c can be understood in terms of recombination. Fragmentation is ruled out. Although 
the calculated results on A/K and Q/4> ratios do not agree perfectly with the data on the 
high end of the pt range, they represent only less than 30% discrepencies in the spectra 
themselves that vary over 5 orders of magnitude. The most significant revelation is that a 
robust thermal parton source is created at RHIC before hadronization. 
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Figure Captions 

Fig. 1. Transverse momentum distribution of in central Au+Au collisions. Data are 
from |23j. The solid line is the sum of the three contributions: TT (dashed line), TS 
(dashed-dot line), SS (dotted line). 

Fig. 2. Transverse momentum distribution of A in central Au+Au collisions. Data are from 
[23] . The heavy solid line is the sum of the four contributions: TTT (dashed line), 
TTS (dashed-dotted line), TSS (dotted line), SSS (light solid line), (a) no valence 
quarks are included; (b) valence quarks are included. 

Fig. 3. The ratio of A to K®. Data are from [221 [23]. The solid line is from Fig. 1 and Fig. 
2(a). 

Fig. 4. Transverse momentum distribution of <fi in central Au+Au collisions. Data are from 
[T9| 120] . The solid line is the sum of the three contributions: TT (dashed line), TS 
(dashed-dotted line), SS (dotted line). 

Fig. 5. Transverse momentum distribution of Q in central Au+Au collisions. Data are from 
[21| [241 [25] . The heavy solid line is the sum of the four contributions: TTT (dashed 
line), TTS (dashed-dotted line), TSS (light solid line), SSS (dotted line). 

Fig. 6. The ratio of Q to <ft. Data are from [211 [25]. The dashed line is for c = 7 = 5 in 
the RF and for thermal parton recombination only. The solid line is for c = 7 = 5 
and includes all components of recombination. The dashed-dotted line is for c = 7 = 4 
with all terms included. 
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